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Abstract 

Signal of lighter top squark has been looked for using Tevatron data in the di- 
electron plus di-jet channel. We find that the mass of the lighter top squark when 
it decays dominantly in the electron plus jet channel, can be ruled out up to 220 
GeV at 95% C.L. using di-electron data. In the framework of R-parity breaking 
SUSY model we exclude relevant R-parity violating couplings for a range of top 
squark masses and other SUSY parameters. The bounds on R-parity violating 
couplings are very stringent for the parameter space where lighter top squark 
turns out to be the next to lightest supersymmetric particle. 



PACS no: 11.30.pb, 14.60. Cd, 14.80.Ly 



I. Introduction 



The Minimal Supersymmetric Standard model(MSSM) [1] so far is one of the most credible 
candidate for the beyond standard model(SM) physics. There is no single evidence of su- 
persymmetric (SUSY) particles, however, no observation can rule it out either. Therefore, 
hunting for SUSY in the next generation of colliders at Fermilab and at LHC experiments 
is a very challenging programme. At present, from the non observation of SUSY signals in 
the past experiments, mainly at LEP [2] and Tevatron [3] , masses of SUSY particles have 
been constrained. 

In MSSM, there is a mixing between the scalar superpartners of the two chirality states 
of fermions, fi and Jr. The extent of mixing of these two chiral states is controlled by the off 
diagonal term rrif^Af — (xtanfl) in the sfermion mass matrix. It is obvious that the sfermions 
which are superpartners of the massive fermions will have larger mixing effect because of the 
explicit dependence on the corresponding fermion mass m/ i.e. sfermions of third generation 
receive a large splitting between two mass eigen states. Thus the two chirality states of top 
squarks £l,£r, has large mixing,resulting in large splitting between the two physical mass 
states ti,t 2 (assume ~ m t ~ 2 ) [4]. Moreover, because of the large Yukawa coupling, the 
soft SUSY masses (m^ , rrif ) also receive a large correction via the renormalisation group 
equation [5] which can push the lighter mass eigenstates, ti, even below the top quark mass. 
Consequently, in a certain region of SUSY parameter space it may turn out to be the next 
to lightest SUSY particle(NLSP), the lightest neutralino Xi being the lightest SUSY particle 
(LSP). It is to be noted that in the canonical SUSY searches at colliders the missing energy 
due to the presence of xi which is assumed to be stable and non interacting, plays a very 
crucial role [5]. 

In hadron colliders t\ can be produced copiously, since it is colored and comparatively 
lighter than the other sparticles. Therefore, in the context of SUSY searches at hadron 
colliders, top squark searches has received a special attention. The search strategy of top 
squark depend very crucialy on its decay pattern. As for example, the loop induced flavour 
changing neutral current decay mode [6], 



yields acoplanar jets and missing energy from top squark pair production. At Tevatron, 
data corresponding to RUN-I experiment has been analysed to find top squark signal in 
this channel. Negative results have constrained lighter top squark mass, m~ tx £119 GeV (102 
GeV) for m^o =40(50) GeV [7]. However, if kinematically accessible, the top squark decays 
dominantly into a lighter chargino(m-±) and b quark, 



Because of the cascade decays of Xi m t° neutralino and massless fermions, Xi ~^ Xiff, the 
top squark pair production leads final states consisting leptons and jets accompanying by 
missing transverse energy [8]. Beside these popular decay modes eq.l and 2, there are also 
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other interesting decay channels which yield a variety of signals in the colliders. All those 
decay modes will be discussed in the next section. 

It is found that top squark mass upto ~ 170 GeV can be probed in RUN-II experiment 
at Tevatron with integrated luminosity 2 fb" 1 per experiments for the entire region of SUSY 
parameter space where t\ state appears to be NLSP [9] leading the decay mode (eq.l) 
with 100% branching ratio. For the scenario, when is heavier than lighter chargino 
mass, 4 fb -1 luminosity is required for the same discovery limit in the dilepton plus missing 
energy channel which is heavily contaminated by top backgrounds. The upgraded RUN-II 
experiment which may deliver high luminosity ~ 20 fb^ 1 may extend this reach upto ~ 220 
GeV [9]. 

However, top squark phenomenology shows up new features in the framework of R-parity 
violating(RPV) SUSY models. In the SUSY models R-parity conservation(RPC) is assumed 
to forbid the decay of proton ensuring conservation of lepton and baryon numbers, L and 
B respectively. However, one can avoid proton decay problem by invoking either L or B 
conservation. Thus one can have two kind of RPV SUSY models corresponding to L or B 
violation. As a consequence, in any SUSY cascade decay process within the framework of 
RPV SUSY model, the LSP can have decay modes either in the leptonic or hadronic channels 
leading multileptons and multijets in the final states with or without missing energy. The 
prospects of SUSY searches at Tevatron in the context of R-parity breaking SUSY model 
has been investigated in great detail [10]. 

In RPV SUSY model, the superpotential is 



W W'mssm I U/, v (3) 

where Wmssm is the superpotential containing yukawa type of interactions giving masses to 
the fermions and corresponds to the potential containing terms which violate L and B 
numbers, 

W % = \ ijk L t LjE K + X'^LiQjDk + X'^J) I), I ),, + mUE^ (4) 

Here i,j, k are the generation indices, A, A' and A" are the dimensionless yukawa couplings. 
The superfields L,Q represent the SU(2) doublets for leptons and quarks respectively where 
as singlets U,D and E stand for up type, down type quarks and charged leptons respectively. 
Last term in eq. 4 mixes the mass terms of the lepton and higgs doublets [1 1] . However, in 
the present case we will work in the context of spontaneous [12] RPV neglecting this bi-linear 
term in eq. 4. 

In this work we focussed on the top squark decays in RPV SUSY model. Considering 
only lepton number violation the X' i3 j coupling which leads to a new decay channel of t±, 

h^l + q- (5) 

As a consequence, in this scenario the pair of top squark production are signalled by dilepton 
plus di-jets. As we know, the identical final states also appear due to the pair production 
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of Leptoquark(LQ) assumed to be the composite object of lepton and hadron [13] and its 
subsequent decays into the corresponding lepton and quark 

LQ-*£ + q. (6) 

At Tevatron experiments, the searches for all three generations of LQ which are signalled 
by ee + 2jets, fifi+2 jets and rr+ 2 jets respectively were carried out [14, 15, 16]. This 
remarkable similarity between the final states due to the pair production of top squarks in 
RPV SUSY scenario and Leptoquark pair production and its subsequent decays via eq.6 
motivated us to exploit the existing data corresponding to this Leptoquark searches to con- 
strain top squark mass which is involved in the production mechanism. Recall that t\ states 
has also other RPC decay modes as well which will be discussed in the next section. The 
branching ratio(BR) of t\ in the RPV channel, eq.5, is controlled by RPV couplings and 
as well as other SUSY parameters which are involved in determining the decay rates corre- 
sponding to RPC decay channels. In this work our main goal is to analyse top squark signal 
in the di-electron plus di-jet channel and compare it with the data corresponding to this 
final state which was used for 1st the generation of Leptoquark search at Tevatron. This 
leads to lower limit of top squark mass for a given BR of top squark in the RPV channel. 
Moreover, the RPV couplings (A' 13j ) can be excluded for a given top squark mass and SUSY 
parameter space. In the case of dimuon channel, we found in the paper of Ref [16] that 
the background corresponding to this channel has been analysed using neural network (NN) 
analysis and the same NN has been trained to analyse signal. Because of this we could not 
perform the analysis for this channel, while the investigation of the ditau channel due to 
the A 33? - RPV coupling is now under progress. Hereafter, whenever we will refer to data, 
that will correspond to the data in the di-electron plus di-jet channel. Earlier there was also 
attempt to constrain RPV couplings using top quark data [18]. 

It may be recalled that these RPV decay channels drawn a lot of attention for the possible 
interpretation to explain the excess of high Q 2 events reported in HI and ZEUS experiments 
few years back [19]. In that context, implications of these channels was examined at Tevatron 
[20]. 

It is worth mentioning here the existing bounds of the RPV couplings which are relevant 
for the present purpose. From the direct production of t x in electron-proton collision at 
HERA, the HI experiment predicts a bound, A' 131 £ 0.05(.02) at 95% C.L for m h =200(100) 
GeV [21]. A bound on A' 131 also exist from atomic parity violation(APV), A' 131 -0.07 at 95% 
C.L for 771,^=200 GeV [22]. The forward-backward asymmetry in e + e~ collision predicts the 
bound on A' 132 ~ 0.28 for m~ tL =100 GeV [23]. The stringent bound on A' 133 comes from 
neutrino data. In the framework of RPV SUSY model the neutrino masses can be generated 
from the tree level contributions due to the bi-linear terms or loop contributions from the 
trilinear A and A' couplings [24] . The detailed phenomenological analysis has been done using 
neutrino data including trilinear and bi-linear couplings [25, 26]. The most favoured large 
mixing angle solution constraints the trilinear couplings A' 133 ~ 10~ 4 assuming M susy = 100 
GeV [26]. 
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We have organised our paper as follows. In section II, we discuss branching ratios of 
top squark decay into various channels in the context of RPV SUSY model. Our analysis 
is described in section III followed by our results in section IV. Finally we conclude with a 
summary in section V. 

II. Top squark decay in R-parity violating SUSY model 

As mentioned in the last section that the lighter state of top squark, t\, has many phe- 
nomenologically interesting decay modes depending on its mass. The most dominant decay 
mode of t±, if kinematically accessible, is via lighter chargino state, Xi an d b quark, (eq.2). 
In the absence of this two body charged current decay mode, the flavour changing neutral 
current decay mode and the four body decay channel, into a b quark, the LSP(x?) and two 
approximately massless fermions [27, 28], 

h - bxlff (7) 

are the only available decay modes and very competitive to each other. This 4-body decay 
channel occurs via many diagrams involving a variety of heavier SUSY particles in the 
intermediate state. In the papers of ref. [27, 28], the decay pattern in this channel has 
been discussed elaborately over a wide range of SUSY parameter space assuming R-parity 
conservation. Surprisingly, in certain region of SUSY parameter space the 4-body decay 
mode takes over the loop level decay mode and the branching ratio may shoot up to ~100%. 
It implies that the signal corresponding to the neutral current decay modes of top squarks 
will be suppressed. The impact of this 4 body decay modes in the context of top squark 
searches at upgraded Tevatron has been discussed [29, 30]. Beside these decay channels of t±, 
there are also a few other decay modes which may be interesting from the phenomenological 
point of view. As for example, if kinematically accessible, the three body decay mode to 
bottom quark and a W boson or a charged Higgs scalar H ± , and a neutralinos, t\ — > bW ± X\ 
or ti — > bH + Xi ma Y open up [31]. Moreover, in the light slepton scenario which is viable in 
some SUGRA models [32] t\ decays via the final states containing sleptons, ti — > btv, b£u [33]. 
In RPV SUSY model, top squark decay channel, eq. 5 opens up due to the interaction 

Xi3j£ihqj + h.c (8) 

which is a subset of the Lagrangian given by eq.4. The species of the lepton and quark 
depend on the choice of % and j respectively, where i, j— 1,2,3. Since we are not restricted to 
the jet flavour therefore j can be of anything 1-3 in the decay process,eq. 5. 

Neglecting the fermion masses, the decay width of ti in the RPV channel (eq. 5) is given 

by 

I>(fi -> I + q) = J — m h (9) 
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9i is the mixing angle in the top squark sector; it appears due to the replacements of ti by 
the physical states t\. The BR in this channel is given by 



e = 



ryti -> i + g) 



(10) 



- i + g) + r& - #pc) 



where r(ti — > RPC) stands collectively for the total decay width in all accessible RPC de- 
cay modes of i\ states. The rate of those RPC decay modes depend on SUSY parameters, 
particularly on the SU(2) gaugino mass M2 (assuming electroweak gaugino mass relation, 
Mi ~ where M 1 is the U(l) gaugino mass), the higgsino mass parameter ji and tan/5 
- the ratio of the vacuum expectation values of two higgs doublets which give up type and 
down type quark masses. Beside these parameters, it involves other sparticle masses and 
parameters, e.g. mass of sleptons (rrig), mass of squarks (rrig) and A-terms, the trilinear cou- 
plings. We investigate the relative rates of the RPV decay mode over the RPC decay modes 
for a wide range of SUSY parameter space. We compute the decay widths corresponding 
to the loop induced flavour changing decay mode (eq. 1) and four body decay width(eq. 7) 
following the expressions given in the paper of Ref. [6, 27, 28]. 

In fig.l, we present contour plots for the constant BR (e ) corresponding to the RPV 
decay modes of ti in the - X' i3 j plane for a fixed set of SUSY parameters: M 2 = 125 GeV 
H = 400 GeV, tan = 4, m q = 300 GeV and m i = 200 GeV and cos 6~ t = 0.8. We find for low 
m^, where the two body chargino decay mode (eq.2) is kinematically inaccessible the RPV 
decay modes has comparatively appreciable rates even for very small value of A^(~ 10~ 4 ). 
In contrast, for higher t\ masses, the two body charged current decay mode will open up 
resulting in the supression of the RPV decay. In this case the RPV mode will be important 



At Tevatron top squark pairs are produced via quark-antiquark annihilation and gluon gluon 
fusion, 



Since, this pair production mechanism is dominated by the QCD process, the cross section 
depends solely on m~ tl [34]. The SUSY-QCD corrections enhance the cross section by another 
~ 30% over most of the SUSY parameter space [35]. We estimate the cross section setting 
renormalisation and factorisation scale at Q 2 = s and use CTEQ3L [36] for the parton 
distribution functions. The typical top squark pair production cross section ranges from ~ 
10.-0.1 pb for m h -100-200 GeV at ^Ts = 1.8 TeV. 

We analyse the top squark pair production in the di-electron plus di-jet channel which 
originates because of the top squark decay via eq. 5. As we discussed earlier that the same 
type of event topology also appears due to the pair production of 1st generation of Leptoquark 




O(0.1). 



III. Squark production: Event Analysis 



QQ, 99 



(11) 
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and its subsequent decay through eq. 6. The dominant SM backgrounds corresponding to 
this signal come from Drell-Yan process. The other sources of backgrounds are due to the 
production of WW, VF+jets, tt and Z+jets followed by the leptonic decay of vector boson. 
For the present purpose we do not estimate the cross sections for each of these background 
processes. Instead we closely follow the analysis as described in the paper [15] for 1st 
generation of Leptoquark searches at DO experiment, where the kinematic cuts are chosen 
to effectively suppress these backgrounds. 

The number of signal events for a given is given by, 

n sig = a hi ,.e d .C.e (12) 

Here, stands for the detection efficiency which includes the acceptance efficiency due to the 
selection and as well as background rejection cuts and all systematic efficiencies like trigger 
and lepton identification efficiencies. The luminosity for the given set of data is given by C 
and cr^^ represents the top squark pair production cross section and e is the BR of t\ decay 
into lepton plus jet. We generate events using PYTHIA(V6.206) [37] producing top squark 
pair which are forced to decay via eq. 5 in order to assume the 100% branching ratio of top 
squark in this RPV channel. Our intention is to compute the accepted efficiency for a given 
set of kinematic cuts corresponding to this di-lepton plus di-jet final states. We take into 
account the effect of initial and final state radiation as well as fragmentation effects in the 
event generation. We adopt the following strategy to compute the signal cross section for a 
given and then comparing our results with the existing data we obtain limits of and 
RPV couplings. 

• First, for a given m^, we compute acceptance efficiency for the signal process by 
generating events using PYTHIA [37] applying the same set of cuts which are used in the 
analysis for the first generation of Leptoquark searches [15]. Then we multiply the respective 
systematic efficiencies, e.g. trigger and lepton identification efficiencies appropriately with 
the acceptance efficiency to obtain overall detection efficiency e^. 

• Secondly, from the existing data corresponding to this final state, for a given e and 
m^, we estimate the cross section limit using eq.12 following Bayesian approach with a flat 
prior probability distribution of cross section. The systematic and statistical uncertainties 
are included taking Gaussian prior distribution of each of them. Then the limit of cross 
section is compared with the theoretical prediction for a given setting a fixed value of e. 

• And finally, repeating this procedure for various m^, the upper limit of e can be 
obtained for each m^. Eventually, this upper limit of e can be translated to obtain upper 
bound of RPV coupling for a given set of SUSY parameters. 

Now, in the following we discuss event analysis for the di-electron plus di-jet channel. 
The pair production of top squark and subsequent decays of each top squarks in the RPV 
channel via coupling X' 13 j 

titl — > ee + qq (13) 

results the final states, e e + 2 jets. As discussed earlier that it is identical to the final states 
due to the pair production of 1st generation of LQ and its subsequent decay via eq. 6. For 
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our signal Monte Carlo, we followed the analysis very closely as described in Ref. [15] where 
search for 1st generation LQ with same final states has been described and then compare 
our results with the data. As we mentioned above that the dominant SM backgrounds 
corresponding to this di-electron final state are from Drell-Yan production with two or more 
jets, ti production and multijet events in which two jets are miss identified as electrons 1 . In 
our simulation we selected events in the hadronic and electromagnetic calorimeter cells in 
pseudorapidity and azimuthal angle(^) of size Arj x A<p = 0.1 x 0.1. Cells with E T > 1 GeV 
are taken as initial seeds to form calorimetric tower. For the jet reconstruction we use the 
routine PYCELL in PYTHIA [37]. Jets are reconstructed with cone radius 0.7 and accepted 
only those which has transverse energy Et > 8 GeV and are smeared by 0.5 x y/Ex- Events 
are subjected to the following sets of cuts mentioned in the paper of Ref. [15] , 

1. Two electrons with E T > 20 GeV and within the coverage for central calorimeter (CC) 
\rj\ <1.1 and for endcap calorimeter (EC) 1.5< \rj\ <2.5. 

2. At least two jets having E T > 15 GeV and \rj\ <2.5. 

3. Isolation between electrons and jets are maintained by requiring AR e j > 0.7 where 
AR = + Ar]' 2 

4. Events of having di-electron invariant mass between 82 < M ee < 100 GeV are rejected. 

5. The total visible transverse energy (Sy) satisfy the cut St > 350 GeV where St = 
H T + H Tl H T = sum of the E T of the two electrons; H T =sum of the E T of all jets. 

The cuts 1-3 are the event selection cuts where as cut 4-5 are the background rejection 
cuts. The cut 4 is used to avoid the contamination due to the events from Z production. 
In the Drell Yan process, leptons and jets are not so hard as in the case of signal process 
for mi x ~100-200 GeV. So a cut on the sum of the transverse energies of visible particles in 
the final state drastically reduce this background. The cut 5 serves that purpose and brings 
down the level of background to a negligible level. We found that the signal acceptance 
efficiencies which is only due to the kinematic cut effect vary from 0.1 - 15% for = 
80-200 GeV. We take into account the electron identification efficiencies which are 74 ± 3%, 
66 ±4% and 68 ±9% for CC-CC, CC-EC and EC-EC regions respectively [15] by multiplying 
appropriately with the acceptance efficiencies which are obtained from PYTHIA, and we refer 
this as a detection efficiency. In table. 1, we show the detection efficiencies folding all other 
systematic efficiencies together for various choices of m^. For a given m^, we can obtain 
from the knowledge of detection efficiency the number of di-electron events for a given e 
and luminosity. The total integrated luminosity is 123 ± 7.0 pb _1 for this di-electron data 
set [15]. 

In the paper of Ref. [15] it is reported that no signal events exist in this di-electron channel 
where as the number of background events after applying all sets of cuts as described above is 

1 The probability of faking a jet as electron is very small ~O(10~ 3 ) [15] 
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0.44 ±0.06 [15]. The uncertainty in background estimation is mainly due to the systematics. 
We exploit this information to obtain the limits of top squark pair production cross sections 
at 95% C.L using Bayesian method for different choices of m~ tx and for a given value of e. 
In this cross section limit calculation we take into account the uncertainties in background 
estimation, in the luminosity measurements and the uncertainty in detection efficiency 

Following the strategy as outlined above, using the data which implies that no signal 
events with a background 0.44 ± 0.06, we obtain limits of top squark pair production cross 
sections at 95% C.L. for two choices of e =0.5 and I. In fig.2 we show these limits (solid lines) 
along with the theoretical predictions(dashed line). In theoretical calculation we multiply 
the K-factor 1.3 with the Born level cross section to take into account the next to leading 
order effect [35]. Notice that, one can rule out top squark mass upto 220(165) GeV for the 
choice of e=l(0.5) in a model independent way. Moreover, as we explained already that the 
upper limits of cross sections which are consistent with di-electron data predict upper limits 
of e for a given m~ tl . In fig. 3, we present these upper limits of e for each value of m~ tl which 
is ruled out by data at 95% C.L. In the context of RPV SUSY model this upper limit of 
e can be translated to obtain an upper limit of respective RPV couplings for a given 
and SUSY parameter space. In fig. 4 we show the excluded region in the A' 13 • — plane 
for a given set of SUSY parameters and for two choices of tan/9 =5 and 30. In each figure 
we excluded region for two extreme values of cos #==0.02 and 0.95. The choice of our SUSY 
parameters for fig. 3 and 4 are (units are in GeV): 

M 2 = 130, /i = 500, tan/3 = 5(30). 
m .± = 515(514), m-± = 121(126), 
m -o = 63(65), rn^o = 122(126) 
m-o = 504(506) m^o = 515(511) 

nig = 300, m l = 200, A h ^ t = 200 (14) 
We discuss results in the next section. 



IV. Results and Discussion 

We have computed the signal cross section in the di-electron plus di-jet channel due to the 
top squark pair production at Tevatron in the framework of RPV SUSY model. In addition 
to the RPC decay modes, t\ also decays via two body decay channel into lepton and quark 
due to the presence of RPV couplings. The relative rates of this RPV decay mode are shown 
in fig.l as contours of fixed value of its BR e in the X' 13j - plane for a given set of SUSY 
parameter space. In the lower region of values, in addition to the RPV decay mode 
(eq. 5) other available RPC decay modes are the loop induced decay channel, (eq. 1) and 
4-body decay mode, (eq. 7) which are of the same order in perturbation theory i.e 0(a 3 ). 
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In this mass region as expected, the RPV decay mode will dominate over the other decay 
modes for most of the parameter space depending on the value of X' 13 j and cos Q~ t . Therefore, 
in this region of m^, even very small value (~ 10~ 4 ) of A' 13? - coupling will yield appreciable 
rates for RPV decay suppressing the two RPC decay modes. Once the value of crosses 
the m~±(=112 GeV) threshold, the two body charged current decay mode (eq.2) opens up, 
which is very much competitive to the RPV decay mode. As a result, large value of RPV 
coupling A' 13 -(~ 0.1) is needed to make the RPV decay mode comparable with the two body 
charged current decay mode for a given m^. We present this result in fig. 1 for a single 
SUSY parameter point. However, we have checked that this pattern more or less exists 
for an entire region of SUSY parameter space which are accessible at Tevatron. We intend 
to emphasise that once the two body charged current decay mode of ti opens up then it 
becomes dominant, otherwise the RPV decay mode is the most dominant one in comparison 
to the loop level and 4-body decay modes. 

In fig. 2 we present the limiting values of signal cross section at 95% C.L for two values 
of e(= 1 and 0.5). We also show the theoretical prediction of top squark pair production 
cross section including K-factor [35] by the dashed line in the same plane. Comparing the 
cross section limits with the theoretical predictions top squark masses can be constrained 
as a function of e. For instance, from the di-electron plus di-jet data we set the limit of 
top squark mass, -220(165) GeV for e=l(0.5). Notice that the limit of depends 
very strongly on e. However, there is a 10-30% theoretical uncertainty due to the choice of 
renormalisation and factorisation scales and parton distribution functions in cross section 
calculations. Note that these limits on m~ tx are obtained in a model independent way. 

In fig. 3, we show the upper limits of top squark decay BR(e) at 95% C.L for various 
values. These upper limits do not depend on any specific models. More precisely, if t\ state 
has the decay channel as eq. 5 then the corresponding BR is restricted by existing data, as 
shown in fig.3. As for example, for 77^ = 100 GeV, the 95% C.L. upper limit on e is 0.35. 

In the framework of RPV SUSY model the decay rate of t\ is mainly controlled by X' l3 j 
for a given and cos%(see eq. 5). Therefore, in this model, the upper limit of e can 
be translated to the upper limit of A' 13? for a fixed and other SUSY parameters which 
determine the decay rates of ti into the RPC decay modes (see eq. 10). In fig.4, at 95% 
C.L, we show the exclusion region in the A' 13? - — m~ tx plane using data. The set of SUSY 
parameters corresponding to this plot is given by eq. 14. Notice that for a lower range of 
m- tx (~ 100 GeV) where it appears to be NLSP, the RPV couplings are restricted to be 
A' 13j % 10~ 4 (10- 3 ). Note that in this region for lower value of cos%=0.02, the RPV decay 
rate is suppressed (see eq. 5) leading to weaker bounds where as for higher values of cos 0f, 
bounds are turn out to be relatively better. However, the limits also become comparatively 
weaker in the higher side of m~ tl where it is heavier than m -± . This is because, when the two 

body charged current decay mode of t\ (eq. 2) opens up, then it becomes very competitive 
with the RPV decay mode leading lower BR for the RPV channel. Consequently, in this 
region, the BR limit constrains only higher side of A' 13 ■ for a given . 

It is obvious from eq. 5, that the higher values of cos# f ~ will yield more stronger limits on 
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\' 13 j. With the increase of tan/?, the two body loop decay and 4-body decay [27, 28] widths 
enhances as the virtuality between and 2 decreases resulting a suppression of e for a 
given , which eventually leads a less constrained region in the A' 13j - — plane. It is clear 
that for the region of ~ m-± i.e when t\ is not NLSP, the bounds on RPV couplings are 
comparatively weaker. 

V. Summary 

We investigate the di-electron plus di-jet signal due to the top squark pair production at 
Tevatron. Identical final states also appear due to the first generation of Leptoquark pro- 
duction. Exploiting the existing experimental data analysed by DO group at Tevatron in the 
context of Leptoquark searches we try to constrain the top squark mass for various values 
of BR of top squarks in the RPV channel. Using DO data and assuming the BR of t\ decay 
via the eq. 5, e=l(0.5) we predict lower limits on -220(165) GeV at 95% C.L. Repeating 
this exercise for various value of we exclude part of the parameter space in the e — 
plane as shown in fig. 3 in a model independent way. 

In the framework of RPV SUSY model, this exclusion region in e — plane converted 
to a corresponding exclusion region in the A' 13j - — plane as shown in fig. 4 for a given set 
of SUSY parameter space. We exclude X' 13 j ~ 10~ 4 for =100 GeV and tan/3 = 5 where 
as for high tan/3=30 region this limit turns out to be relatively weak, X' 13j Z 10~ 3 . Notice 
that when t\ state appears to be NLSP, the limits are very stringent and comparable to 
the limit, for the case A' 131 , obtained from neutrino data [26] as discussed in Sec. I. However, 
for, mi x ~ the present analysis does not give any better limit than the others obtained 

from HI experiments and APV measurements [21, 22] and also from neutrino data [26]. We 
conclude that our predicted bounds are very stringent in the region where t\ state appears 
to be NLSP. 
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m fl (GeV) 


Detection efficiency(%) 

e+e- 


100 


2.0 


120 


4.6 


140 


9.2 


160 


15.0 


180 


21.8 


200 


26.7 


220 


29.8 



Table 1: Di-electron plus di-jet detection efficiencies for various m~ h 




Figure 1: Branching ratio contours for the decay channel t 1 — > £ + q. The SUSY parameters 
are: M 2 = 125 GeV, n = 400 GeV, tan (3 =4, cos0 t ~ =0.8 and m~ q = 300 GeV, m i = 200 
GeV, A b = A T = 200 GeV. 
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Figure 2: The top squark pair production cross section limits at 95% C.L (solid lines) 
for e—1 and 0.5 along with the theoretical prediction(dashed line). 




Figure 3: The excluded region by di-electron data at 95% C.L. 
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Figure 4: The excluded region(hatched) by di-electron data at 95% C.L. The SUSY param- 
eters are: M 2 =130 GeV, /i=500 GeV, m- q =300 GeV, m| =200 and A-terms=200 GeV. 
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